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ABSTRACT 


The nonlinear coalescence instability of current carrying 
solar loops can explain many of the characteristics of the solar 
flares such as their impulsive nature, heating and high energy 
particle acceleration, amplitude oscillations of electromagnetic 
emission as well as the characteristics of 2“D microwave images 
obtained during a solar flare. The physical characteristics of 
the explosive coalescence of currents are presented in detail 
through computer simulation and theory. 

Canonical characteristics of the explosive coalescence are: 
Cl) a large amount of impulsive increase of kinetic energies of 
electrons and ions, 

(2) simultaneous heating and acceleration of electrons and ions 
in high and low energy spectra, 

(3) ensuing quas i-pe r i od i c amplitude oscillations in fields and 
particle quantities, 

<4) the double peak (or triple pe ak) s t rue tu r e in these profiles, 
A single pair of currents as well as multiple currents may 
participate in the coalescence process, yielding varieties of 
phenomena. 

In particular, double sub-peak structures in the qu a s i -p e r i o d i c 
osillations found in the time profiles of two solar flares on 
June 7,1980 and November 26, 1 982 are well explained in terms of 

the coalescence instability of two current loops. This interpre- 
tation is supported by the observations of two microwave sources 
and their interaction for the November 26,1982 flare. 
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1. Introduction 

The solar flare phenomenon is a manifestation of an explosive 
release process of energy stored in the lower co rona, i nvo 1 v i ng 
the plasma heating up to 5xlO^K,the acceleration of charged 
particles up to the order of the rest mass energy of electrons 
and ions, as well as the production of electromagnetic radiation 
from the radio band to T -ray wavelengths, (for a previous summary 
of solar flares, see Svestka, 197G;Sturrock, ed. , 1980;Priest, 1982) . 

After launching of the Solar Maximum Mission (SMM) and 
Hinotori satellites, it becomes clear from the observations of 
hard X-rays and T-rays that the electrons with energies up to MeV 


as ions with energie: 


IP to GeV can 


s imul taneous ly 


accelerated within a second during the impulsive phase of a solar 
flare. In a particular flare (03: 12UT of June 7,1980)the T -rays 
showed a quas i — pe r i od i c amplitude oscillation which 

closely correlated with the quas i— pe r i od i c i ty in both microwave 
bursts and hard X-ray bursts. 

Direct observations in soft X-rays (Howard and Svestka, 1977) 
showed that in the active regions there exist multiple coronal 
loops which might carry plasma currents. The interconnecting 
coronal loop might be a quite important physical process for 
energy release in the solar corona. 

Recently, the observations of the interacting coronal loops 
which lead to solar flares have increased in various kinds of 
observations from Hof(Kurokawa et al.,1985), radio (Nakajima et 
a 1 . , 1 984, Kundu, 1985) as well as hard X“ ray (Machado, 1985) . 

In order to explain the rapid quas i— pe r i od i c particle 

acceleration of both electrons and ions observed in the June 

7,1980 flare, (Taj ima et al. , 1982, 1985) showed that the most 

likely mechanism for the impulsive energy release in solar 

flares is the current loop coalescence instability (Finn and 
Kaw, 1 977, Pr i tche t t and Wu,1979). 

It has been shown that by simulation and theo ry (Ta j ima and 
Sakai, 1985) that during the coalescence of two current loops, the 
magnetic energy stored by the plasma current can be explosively 
transformed to the plasma heating as well as the production of 
high energy particles within a transit Alfven time across the 
current loop (which is about 1^10 seconds for appropriate radius 
of the loop) through the magnetic reconnection process. 
Furthermore, the energy release can be achieved with quasi- 
periodic fashion whose periodicity depends on plasma parameters 
such as plasma beta ratio (^), the ratio Bp/B-r between the 
magnetic field. Bp produced by the plasma current and Bt the 
potential magnetic field, as well as the colliding velocity of 
two current loops that is determined mostly from the initial 
total plasma loop current. 

The plasma can be heated up to 60 times of the initial 
temperature and the electrons and ions can be accelerated 
simultaneously by the transverse electrostatic field which can be 
produced during the explosive coalescence process. 

The present paper is to show results obtained from current 
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loop coalescence plasma dynamics with applications to solar 
flares. The current loop coalescence itself is similar with Gold 
and Hoyle model (1960) and emerging flux model (Hayverts et al. , 
1977) in its morphological spirit of the model. However, the 
details are quite different from each other. In particular, 
Taj ima and Sakai (1985) presented the physical basis of such a 
morphological model, i.e., the physics of magnetic reconnection 
process, as well as particle acceleration mechanism. We found 
basically two types of magnetic reconnection, namely, slow and 
explosive processes: the reconnected magnetic flux d'}' is 
proportional to t^(m>l) the slow reconnection process and 4*1^ is 
proportional to (t<,~t)”"' (m>0) , in the explosive process. We also 
found in the ensuing stage the qua s i-pe r i o d i c reconnection with 
qua s i “pe r i 0 d i c acceleration during the coalescence of current 
loops. 


2. Physical Characteristics of Current Loop Coalescence 

In the problem of solar flares, the reconnection of field- 
lines is believed to take place due to finite resistivity. The 
paper based on a boundary layer analysis by Furth, Killeen, and 
Rosenbluth found an instability (the tearing instability) which 
grows at a growth rate T =C Ta 

Sweet and Parker, on the other hand, obtained a steady-state 
solution which has a narrow x-point angle; the time scale T$p of 
reconnection is characterized by "Tsp = ^ • whe r e 
2L is the length of the impinged plasma, and are the 
densities inside and outside of the singular layer. Petschek 
similarly obtained a steady-state solution which has a large x- 
point angle; the reconnection time scale Tj, is independent of 
res i s t i V i ty,'lj,"'ot .The time scale of reconnection due to^ the 
mechanism of the Sweet-Parker process is characterized by6‘-*C6 = 

, wh i 1 e that of Petschek is by 6^ . It may be said that the 
paper by Furth et al. is appropriate for situations of 
spontaneously growing tearing modes, while the papers by Sweet et 
al. are for problems of driven reconnection, although the 
distinction is yet to be clarified in more precise scientific 
terms. 

Since the paper Furth et al. was published, a considerable 
amount of literatures have been written which further 
investigate the nonlinear processes of the tearing instability. A 
paper by Rutherford discussed nonlinear secular growth of a 
single tearing mode, while a paper by Drake and Lee found the 
collisionless equivalent to the Ru t he r f o r d^ r e g i me . The Rutherford 

time scale Tft is characterized by T;j— where C is a 

constant less than unity, is the shear field, and 'f' the 
reconnecting flux. Carreras, Rosenbluth, Diamond et al. discussed 
nonlinear stages of many tearing modes with different helicities 
(i.e., on different rational surfaces), emphasizing mode 
couplings. With more than one helicity, the coupling is 
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inherently three dimensional. 

It has been recognized that the processes of magnetic 
reconnection are quite different depending upon whether 
reconnection is driven or spontaneous. For the spontaneous cases, 
as we briefly reviewed in the above paragraph, the current sheet 
becomes unstable against the tearing instability evolving from 
the unstable magnetic configuration. It is customary to argue 
that the Rutherford regime follows the linear stage and then the 
highly nonlinear mode coupling stages come into play. There are 
natural disturbances such as solar flares and geomagnetic storms, 
in which some of the observed time scales of the entire process 
of explosive phenomena are as small as the Alfven time. The 
scenario for reconnection based on spontaneous tearing has to 
confront the observed rapid reconnection processes, although the 
nonlinear stages of spontaneous reconnection can be quite rapid, 
as reviewed in the preceding paragraph. In this case, one has 
to maintain the unstable configuration in equilibrium long enough 
to reach the disruptive stage. On the other hand, driven 
reconnection is much faster and therefore quickly enters the 
nonlinear stage. This is one of the reasons why the driven 
reconnection has been studied by many authors. 

In order to rapidly convert magnetic energy into kinetic one 
by a substantial amount, it seems necessary that the bulk of 
magnetic energy has to participate in the conversion process: 
the resistive heating at the x-point alone is to meager. This is 
because the available magnetic energy at the X“point is small by 
itself. On the other hand, the ideal MHD instabilities such as 
the kink instability and the coalescence instability are the 
processes that involve the bulk current redistribution in a 
matter of the Alfven time scale. 

In the present paper we pick the coalescence instability as 
the primary instability to investigate its nonlinear 
consequences. we do so because (1) although it is an ideal MHD 
instability for drive in the linear ser\se, it would not 
nonlinearly evolve if there is no resistive Cnon-ideal MHD) 
effect; (2) therefore, it can involve a large amount of 
conversion of magnetic to kinetic energies in a short time; (3) 
it is essentially a two-dimensional instability, thus more 
amenable to thorough analysis of the fundamental processes of the 
instability. It is interesting to observe that with this two- 
dimensionality restriction we still find an explosive process as 
we shall see. The coalescence instability starts from the Faddeev 
equilibrium, which is characterized by the current localization 
parameter 6c • The equilibrium toroidal current (in the z- 
direction) is given as =B^^k ( 1 - 6^) (c o sh ky + c o s kxT .The para- 
meter varies from 0 to 1 with small corresponding to a weak 
localization and close to unity corresponding to a peaked 
localization; in the limit of^c-^1 the current distribution 
becomes delta f unc t i o n- 1 i ke . The rate of reconnection was that of 
Sweet-Parker for small<^c» while the reconnection rate experiences 
two phases for 1 a r g e r ^bu t smaller than 0.8). This emergence of 
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two phases is similar to the case of the driven reconnection. The 
intensity of coalescence and the rate of subsequent reconnection 
are controlled by just one parameter, the current localization 
C 6c ) • In this problem there is no ambiguity as to the nature of 
the driver in contrast to the reconnection driven by external 
boundary conditions. For the case with =0.7, the second phase 
showed the reconnected flux increasing as t with l<c^<2. 
(Brune l“Ta j ima scaling) This indicates that the more the current 
localizes, the faster the reconnection becomes. 

As seen in the next section, for the case with 6^=0. 85 we found 
more rapid reconnection called explosive reconnection, in which 
the reconnected flux/i^ increases as C t^ - t)'*^'^^ . (Ta j i raa“Saka i 
s c a 1 i ng) 

2. 1 Summary of Simulation Results 

We combine both a MHD model and a kinetic model of 
simulation. The results from these two different models are 
consistent in basic points, but are complementary in many 
detailed aspects. The kinetic simulation model we adopt here is 
the electromagnetic particle code with 2“ dimensions. 

^ The MHD simulation model we use is the MHD particle code with 
2-2 dimensions. The MHD particle code is robust in applications to 
problems even with strong turbulence, flows, convections, and 
density depression. This is helpful because the present problem 
involves fast (explosive) reconnection, strong density depression 
and compression, and strong flows. 

Results from the electromagnetic particle model are now 
discussed. Figure 1 displays the typical time history of various 
field and particle quantities observed in our simulation in which 
after the initial transient (up to t=4i^;*in the code unit to be 
explained in the following section) the phase of coalescence of 
two magnetic islands commences. It is seen in Figs. l(a)-(c) that 
around t=27 the magnetic and electrostatic field energies shoot 
up explosively as well as the ion temperature in the direction of 
coalescence (the x“d i r ec t i on) . The unit of computational time is 
omitted hereafter whenever it is unambiguous. It is also seen in 
Figs. l(a)-(c) that (1) after the explosive increase of the field 
energies and temperature this overshooting results in synchronous 
amplitude oscillations of all these quantities with the period 
being approximately the c o mp re s s i o na 1 Alfven period; and (2) 
superimposed on these overall amplitude oscillations is a 
distinct double-peak structure in the electrostatic field energy 
and the ion temperature. Although we are interested in analyzing 
the entire episode of the run including the initial phase and the 
post-explosive phase, we focus particularly on the explosive 
phase of the c o a 1 e s c e nc e . We replot Figs. l(a)-(c) into Figs. 1(d)- 
(f) to find the way in which these quantities increase toward the 
catastrophic point. We find from Figs. l(d)-(f) that (1) the 
magnetic energy explodes as (t^ -t)”^^ ; (2) the electrostatic 

energy explodes as (t^ -t)"^ ; and (3) the ion temperature in the 
coalescing direction explodes as (t^ — t)'^-^^ until saturation due to 
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1. Explosive increase of field energies and temperature during the coalescence of two 
magnetic islands: ENf particle simulation results. For the case = 0.2. Other 

parameters are given in Sec. II. Toward the same explosion time t = 27(i7“*), 
the magnetic energy J5“ (a), electrostatic energy El (b), and the ion temperature 
in the i-direction Tix (c) diverge as (/„ - - 0”^* as 

shown in (d), (e), and (f), respectively. (\Vc took ^ 1.6 x 10^, -- 7.5 x lO’, 

and Tizo ~ 0-85 for the pre-explosive phase values; sec Figs. l(a)-(c)j. 


398 







2. Temporal profiles of particle and field quantities for the coalescence process ob- 
tained from the EM particle simulation. = 0.2. (a) The thick line represents 
the magnetic energy, the thin one the ion temperature in the i-direction. (7“,i at 
t = 0 was 0.5). (b) Electrostatic field energy in time, (c) Ion temperature in the 
2 -direction, (d) Electron temperature in the x-direction. (e) Electron temperature 
in the 2 -direction (f) Inductive electric field (E^) in time. 
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overshooting sets in, where t^, is the explosion time measured 
here to be t^-^27 C d};') in this run. Figure 2 presents the time 
history of various field quantities and temperatures in the 
course of the early formation and the coalescence with the 
toroidal field being such that tBz /m^ c=5tt=0. 2o^. In Fig. 2(a) both 
the magnetic field energy and the ion temperature in the 
direction of coalescence Cx) show that after the early (t<3 ) 

rise which corresponds to the magnetic islands formation by 
tearing a long relatively dormant period (t=3~22) sets in, 
followed by a stage Ct=22~27) rapid and huge increase in these 
quantities. It is also evident that after the rapid increase 
Ct>27) salient amplitudR oscillations ensue due to overshooting. 
It is to be remarked that all the other quantities shown in Fig. 2 
(a)-Ce) closely follow the pattern of Fig. 2(a) with their 
characteristic events simultaneously occurring. It is also noted 
that the amplitude oscillations of the temperatures (Tc* and T;zas 
well as Ti* and T«z) and the electrostatic field energy have a 
structure of marked double peaks. The valley of the double-peak 
structure coincides with the peak of the magnetic field energy 
amplitude. As mentioned in Sec. 1, it is important to notice that 
the rapidness of the increase of each quantity differs and that 
each quantity explosively increases characterized by a certain 
definite, but different, index of exp 1 o s i on (i . e . , t he exponent to 
the time measured backward from the point of explosion time) 
until the saturation stage sets in. The stop of rise (t 3) of 
each quantity in Fig.2(a)-(e) corresponds to the completion of 
islands formation. The following quiescent period (3<t<20) 
corresponds to the stage where the formed islands slowly attract 
each other. The rapid explosive rise (t>20) marks the commence- 
ment of the explosive coalescence. The following stage of ampli- 
tude oscillations correspond to the "breathing” (or pulsations) 
of coalesced islands (c omp r e s s i ona 1 Alfven oscillations). The 
induced electric field explosively increases when there is 

rapid flux reconnection during the explosive coalescence and then 
oscillates as the magnetic flux in the coalesced island is com- 
pressed and decompressed. 

Some of the above findings can be given by a qualitative 
explanation. Since we shall discuss the explosive phase in 
greater detail later, wc try to pay attention to the amplitude 
oscillation phase in particular here. Once two current blobs 
coalesce, they are bound by common magnetic flux. The larger 
coalesced island continues to vibrate. Within the coalesced 
island the colliding two plasma blobs cause turbulent flows which 
dissipate their energy quickly into heat, thereby reducing the 
amplitude oscillations of temperatures and field. As a result as 
we shall see, the momentum distribution of plasma electrons and 
ions exhibit an intense bulk heating and acceleration of the 
tail. The heating in the poloidal direction (x) is due to 
adiabatic compression and decompression of the coalesced current 
blobs. The eventual bulk heating is a result of turbulent 
dissipation of counte r s t reami ng instabilities. The heating in the 
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toroidal direction is due to heating/acceleration by the 
inductive toroidal electric field which is several times the 
classical Dreicer field and the Vp xB acceleration. The double 
peak in the time development of the temperatures occur just 
before (t = t,'^27) and after Ct = t3''29) the maxima of magnetic field 
(t=t 2 28). In Fig. 3 schematic sequential pictures of plasma 
dynamical behavior during coalescence are shown. At t=ti the 
magnetic (JxB) acceleration of ions becomes maximum so that the 
magnetic flux the behind the colliding plasma blobs as well as 
the plasma blobs are themselves strongly compressed. This plasma 
compression causes the first temperature peak at t = t, . After this 
maximum acceleration phase ions aquire substantial velocities 
along the direction of collision so that they detach from the 
magnetic flux against which ions have been compressed. This 
result in an expansion phase (t=ti) of ions, and hence in an 
adiabatic cooling of the plasma as the magnetic fields obtain 
maximum values. The process reverses after the maximum of the 
magnetic fields at t = t^.^ 29ili* , which gives rise to the second 
peak of the temperature. 

The high energy tail particle acceleration of ions and 
electrons may be qualitatively discussed here. The tail formation 
is probably due to a combination of localized electrostatic field 
acceleration across the poloidal magnetic field and magnetic 
acceleration of the poloidal to toroidal directions. 

Electrons are magnetized and are carried away with the 
accelerating magnetic flux, while bulk ions are accelerated by 
the JxB force. On the other hand, the high energy ions are 
produced and dragged by the charge separation created near the 
compressed flux. The difference of motions between ions and 
electrons around t=ti causes a strong localized shock~like 
electrostatic field, E,, , which propagates with a phase velocity 
of the structure vp =v^ . This Vp xB acceleration causes the 
formation of high energy particle in the toroidal direction. By 
this acceleration process, ions and electrons are accelerated to 
relativistic energies in opposite directions along the toroidal 
magnetic field. 

Next we present the simulation results of the case of multiple 
coalescence process (Fi g. 4) . A number of islands are induced as the 
crowbar current is turned on. An interesting phenomenon here is 
that these islands begin to coalesce one by one to form slightly 
larger ones. In turn, these larger islands coalesce further into 
still larger islands. Finally, we are left with only two (big) 
islands that more resembles two sheets of plasma. Eventually 
these two islands also coalesce into but one. In this multiple 
coalescence process, the early stage presents many less regular 
or spiky time profiles of fields and particle quantities. As time 
goes on, larger and more organized time profiles of these 
quantities come to dominate. Eventually, large scale oscillations 
are observed in Fig. 5 , which are set off after the final two 
islands coalesce into one. This is, therefore, a similar process 
to the previous two islands coalesce process. In this particular 
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choice of parameters heating and acceleration of particles are 
more intense than the previous case. Corresponding energy 
distributions Cspectrum) show a longer and more populous 
ta i 1 . i nd i cat i ng a harder spectrum. However, the basic feature of 
two (for perpendicular) and three (for the z-di rectional 
distribution) distinct distribution regimes is still unchanged, 
implying the same heating/acceleration mechanism at work. When 
we inspect the electron z-direction distribution, the highest 
energy of the tail oscillates as the overall coalesced island 
vibrates. On the other hand, the maximum ion energy of the tail 
in the z-direction, although also oscillating, is less sensitive 
to the island vibration. 

Results from the magnetohydrodynamic particle models are 
presented hence. Fig. 6 shows kinetic energy and the reconnected 
flux upon coalescence as a function of time for the case with 
^=0.85. A theoretical curve (to~t)"^is superimposed on the 
simulation result. During the phase of the rapid increase of 
reconnected flux (t=50-9QA c^') the simulation result matches 
reasonably with the theoretical curve. Beyond t=90 ^ c*' the 
increase begins to be mitigated due to a saturation effect (the 
flux depletion). 

Figure 7 displays the case with ^t=0. 7. The reconnected flux 
increases rapidly wi th A*!* o< t** (m*l. 9) . It is, however, less rapid 
than the case with 6c =0.85. The released energy is also less in 
present case. The case with 6^. =0. 3 was treated, where /i'I'oc t*” with 
•w-1* Thus, it is clear that as increases, the process of 
reconnection becomes faster, changing from the Sweet-Parker rate 
to the Brune 1-Ta j ima rate to the explosive rate. It is also 
remarkable that the explosive increase of reconnected flux during 
the coalescence is observed in the MHD simulation as well as in 
the kinetic simulation discussed earlier. 

The structure and its evolution of the plasma and magnetic 
fields during the coalescence are examined. 

These signatures are consistent with our hypothesis (a) that 
the reconnection takes place by the mechanism of Sweet and 
Parker for coalescence with 6c =0.3. The signatures found in 
6c=0. 3, on the other hand, imply that the reconnection process is 
not that of Sweet and Parker. It shows instead that (1) the 
reconnection angle at the x-point has enlarged ; (2) a 

high density spot near the x-point is formed ; (3) the 

reconnected flux increases faster than the Sweet-parker process 
( A'f'ot t"*” with m 2). These are consistent with our further 
hypotheses (b) that the reconnection is through the process of 
Brunei, Tajima and Dawson for coalescence with 6c =0.7. Later 
(t=160), the system approaches saturation when most of the flux 
available has reconnected as seen in 6c =0.7 (at t = 140 for this 
frame) and the high density region shifted from the x-point to 
the coalesced island hedge. The flow is randomized. 

To supplement our kinetic simulation in order to make 
comparison with the MHD particle simulation, we present a run 
starting from the Faddeev equilibrium. In this case the 
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Temporal profiles of the fluid energy (a) and the reconnected flux (b) for f j = 0.85 

obtained from the MHD particle simulation. A solid line in (b) is a theoretical 
curve discussed in Sec. III. 



80 120 
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Temporal profiles of the fluid energy (a) and the reconnected flux (b) for €e = 0.7 
obtained from the MHD particle simulation. A solid line in (b) is a fit to the 
observation. 
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equilibrium is the same as the MHD particle run, although the 
parameters such as the collisionless skin depth are different 
because of the difference in model. 

Figure 8 presents the pattern of the plasma and fields of the 
case =0.85, where we see faster and explosive reconnection 
(Fig. 6 ). We are advancing our third hypothesis (c) that the 
coalescence with =0.85 is explosive whose reconnection process 
is to be characterized by the present paper. In frames of Figs. 8 
(a) - (d) (t=50) one sees the coalescence behavior before it becomes 
explosive. Although, in Figs. 8(a) and (b) , in particular, one can 
detect some deviation from the Sweet-Parker type for =0.3, it 
is qualitatively similar to the^c =0. 3 case and the =0. 7 case 
at this stage. In Figs. 8 (e) and (f) (t=75) , we now see 
significant deviations in pattern from the cases with less . A 
much wider reconnection angle than the previous ones is observed 
in Fig. 8 (e) . From these observations it can be argued that the 
widening of the reconnection angle has to be accmopanied by fast 
or explosive coalescence. 


2. 2. Theo re t i cal Model of Explosive Coalescence 

As shown in the previous section, the current sheet of nearly 
one-dimensional structure is formed in the explosive stage of the 
coalescence instability. As the coalescence proceeds further, the 
magnetic field structure approaches an x~type (Petsheck type) 
configuration. We assume th a t Vj*. »Vi'i , ViZ , in which x is the 
direction of coalescence, while y is the direction of poloidal 
magnetic field line and z is the direction of plasma current. We 
treat the external plasma dynamics of the explosive stage as a 
one-dimensional problem. 

As we shall see in our separate paper(Tajima and Sakai, 1985), in 
one-dimensional limiting case of two— d imens i ons , we get 
essentially the same results as we obtain in this section. 

We start from the two-fluid model equations of plasma and the 
Maxwell equations, neglecting the displacement current. We assume 
the adiabatic law of states for both electrons and ions. The 
basic equations read as follows: 

^+V-(nyVy)=0 


m 


>n,^=nyey(E+^xB)-Vpy, 
\ ^ 

V X B = — / 
c 

3 

V - E = 47r^njC_,-, 


las 

V X E 57> 

c ot 


(2) 

(3) 

(4) 

(5) 
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0 128 0 X— 120 

8* Spatial structure of plasma and fields “before” and “during” coalescence with 
£c = 0.85. (a) Magnetic field lines, (b) Plasma density contours, (c) Plasma flow 
velocities, (d) Current density (J^) contours, (a)-(d) at t = SOAcJ^, (e) Magnetic 
field lines, (f) Plasma density contours, (e) and (f) at t = S7,SAcJ^. 
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where j denotes the species of particles and is the ratio of 
heat capacity which is related to the degree of freedom of the 
system f as JT =1 +yf. 

During coalescence, there is no specific scaLe length. The 
scale length characterizing the current sheet varies continuously 
in time without deformation of global structure of current sheet. 
If one looks at the evolution of the system locally in time (in 
the neighborhood of t = t^ ) , the system undergoes the rapid field 
and temperature swelling, compression of plasma and the size of 
the current sheet, change of the reconnection angle etc. in a 
certain specific fashion which was detailed in Sec. 2. If one 


looks at the same system locally in time a little later (in the 
neighborhood of t=t^ +t^ ) , the system undergoes these changes with 
different magnitudes, but still in the same specific fashion. 
That is, the relations (laws) that govern the explosive 
coalescence themselves are invariant under the change of time 
scale. This was the manifestation of the presence of 
se 1 f -s imi lar i t y in the system during explosive coalescence. 

Such a physical situation may best be described by self- 
similar solutions in which scale factors vary continuously. We 
introduce scale factors a(t) and b(t) as follows, 

d 

Vex = -X, C7) 

CL 


where a dot represents the time derivative. An ansatz is imposed 
here that the velocities are linear in x. The linear dependence 
on X of the velocity implies that ions and electrons stream in 
the opposite direction around the center of current sheet, x = 0. 
The scale factors a and b will be determined from the above basic 
equations. From the continuity equations of electrons and ions, 
Eq. (1) » we obtain 

n* = no/a, C9) 


n; = rio/6, (10) 

where is a constant. Equations (9) and (10) show that the 

densities of ions and electrons are nearly homogeneous in space 
and vary only in time during coalescence. The self-similar 
solutions obtained here are local solutions in space whose 
physical process dominates near the current sheet. We therefore 
neglect the higher order terms in space proportional to x"* and 
higher hereafter. The current Jz in the sheet is nearly 
constant. This means that as n] is nearly constant, is also 
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approximately constant in space. 
Eq. C3) , we obtain 

Bo(0 _ 4jreno 
"1 c 


(?-?)■ 


Neglecting the term with x* in 
( 0 ) 

\ 

cm 


where we assumed the magnetic-field varies as B^=B<,Ct)^, where 

% is the magnetic field scale length. This ansatz is consistent 

with the assumption that the sheet current is constant in space. 

From the y-component of Eq. C5) and the z-component of equation 
of motion for electrons Eq. (2) we obtain 

C12) 


J3o ” ^ ’ 


i’ O BoW J _ 0 
Xc ^ ■ 




at 


me 


C13) 


C14) 


where 


E. = E,o{t) + 


Equations C12) and (13) yield 

^o(0 = ^, (15) 

c 

where B^o is a constant. From the z— component of equation of 
motion for ions, we get 


dv) 


( 0 ) 


at 


= -^J5e0. 
mi 


From Eqs. (14) and (16) we have 


J») = 


mi 

From Eqs. (15), (17) and (11), we get for v^’ and Ea 


cjBqo^ 


t,(o)= 


E.= 


(16) 

(17) 

(18) 
(19) 


Bopax^ mecBpo d / b ^ 
ea^X 47rno«*A \a*(| + ^) / 

Assuming that the electrostatic field Ex varies like Ex=Eo(t)xA, 
we obtain from Poisson’s equation (4) 


Ep — incnpX 


(i-i)- 


(20) 


Furthermore, the equations of state for electrons and ions give 
rise to 
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(21) 


p _ _ -^o« 

■» e “ 


P = — - 
•» » 


oTf 2a*»+2 A2 ’ 

Pqx •P Ot 


(22) 

5'T 26^+2 A2* 

We now go back to the x-component of equations of motions for 
electrons and ions in order to obtain the basic equations for 
a (t) and b(t). If we neglect the small terms of the order of the 
mass ratio m«./mL , we obtain 

O = -^pe(l ~ ) ” 4nmcnoX'^a‘ men«A2a'» ^23) 


- 2 . Pq{ 

^ ”* a' ^ m,noA26^ 

Furthermore, assuming that the plasma is quas i -neu t ra 1 n =n 
i. e. , a = b, by adding Eqs. (23) and (24) 


(24) 


a = — 


A2a2 A^aT’ 


(25) 


where 

j Bqo j 2 {Poe + Poi) 

Va = r and c* = ■; r — . 

^ 47rno(m,- + rrie) \rn^ + m,)no 

In Eq. (25) the first term of the RHS corresponds to the J x B 
term. This is the term that drives magnetic collapse. The second 
term corresponds to the pressure gradient term. This term may 
eventually be able to balance the magnetic compression (collapse) 
when ^=3. The condition JT = 3 means that the plasma compression 
takes place in a nearly one~d imens i ona 1 fashion so that the 
degree of freedom of the system ^ becomes unity. When iT = 3, we 
obtain 


When = 2(f = 2) . 


on the other hand, 


a = 


- cl) 

A2o2 • 


we obtain 


(27) 


Once the behavior of the scale factor a(t) is determined from 
the above equations, we obtain various kinds of physical 
quantities as follows, in the quas i ~neu t r a 1 plasmas, and 
neglecting the mass ratio (me /mt-^0). 
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( 30 ) 
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no 

m = lie = — 

a 



(33) 

where the electrostatic 

field Ey in the 

quasi ~neu t ra 1 

p 1 asmas 

i s 

determined from the 

equation of motions for ions, 

no t 

from 

Poisson’s equations. 





From Eqs. (28) and 

(29) we find as 

important result that 

i n 


the explosive phase (a 0) the electrostatic field Ex" (a''^+ a"’^) 
grows more rapidly than the magnetic field (B 3 <x a”^) does. This 
effect comes into playing a pivotal role for high energy particle 
produc t ion. 

Now we investigate the global time behavior of coalescence by 
making use of the first integral of Eq. (26). Equation (26) may be 
rewritten as 

« (34) 


where the effective (Sagdeev) potential V Ca) is given 


Dy 




(35) 


where the 
potential” 
schematic graph of the effective 
The value a which satisfies V(a») 


first term may be reminiscent of the ” gravi tat ional 
while the second of the "centrifugal potential.” The 

potential is drawn in Fig. 9 . 
= 0 i s given by 




(36) 


The minimum of the potential, i obtained from ’^V/a a = 0 is 


V • — 

» mm — 


-V’ 


(37) 


2A2cr 

at a = 2a| = c> /v^^ . The potential becomes deeper when the ratio of 
the kinetic to magnetic energy densities (5 = c^ 
zero. This means that 
compared with the pressure 
given by 


/VA approaches 
the driving force J Bis dominant 
term. The first integral of Eq. (34) is 




X^a 




(38) 


where £, is the initial (Sagdeev) "energy” (dimension: 1/time* ) 

in space of stretching factor a. 
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<39) 


t =al- 


2v\ 

UTo 


the 


As seen from Fig. 9 , the explosive magnetic compression 

corresponds that the scale factor aCt) rapidly changes in time by 
orders of magnitude and vanishes. We may call this an explosive 
magnetic collapse. Such an explosive collapse can be realized (1) 
when the effective potential has a sharp and deep potential well 
and this means that ^ /v^ is very small; (2) when the initial 

total energy /2 is nearly zero. On the other hand, whenf/2 is 
close to Yfi>;n t we have oscillations near the potential minimum 
and no explosive collapse. 

If the total energy is given in Eq. C38) , we can find the 
period T^s of nonlinear oscillations by integrating Eq. (38) 

rr - /**’ Oda 

OS 


-r 

•'a,! 


= 27T 




1/2 


where a«i , atx 


and tA 
period as 


£Z/7X2 'A 

are roots of the equation which gives a = 0: 
2 

= 0 , 


2 

a* + ■' ^ a — -■ * 


(40) 


( 41 ) 


In the limit of £ ~ , we find the minimum 


r„,i„ = 211/3^^1^. 


(42) 


Equation (40) indicates that the period T^s of nonlinear 
oscillations becomes longer when £. tends to zero. 

Let us examine the time history of various physical quantities 
based in the qualitative time behavior of a(t) derivable from the 
effective potential V (a) . The magnetic field energy is 
proportional to Bj, which is given by 


B 


2 

y 


— -^00 


If the scale factor a becomes smaller, By must increase. The 
maximum is given by 


dBl 

dt 


= 0 , 


which yields a =0, namely a = a ti • After the maximum, Bj 
decreases again and reaches minimum at a = a-tz • The oscillatory 
behavior of the magnetic field energy is schematically drawn in 
Fig. 10 . The period of the oscillation is given by Eq. (40). 

The electrostatic field E is given by Eq. (29). The time 
history of the electrostatic field energy, which is proportional 
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to , is analyzed by investigating 

dt 

This condition is equivalent to 

^0 = 0 , 


C43) 


0 r 


d£o 

dt 


= 0 , 


where £© (t) = first condition £* = 0 occurs at 

aj = Pe. /min^'vX = cf/v^. The second c ond i t i on 3£e/it=o g i ve s 
conditions, namely 

Ci) a = 0, which occurs at 


(44) 

a = 
two 


a = 


, 3t2 


(i i) 


a = o-i = 


4 Poe ^ 1 ft , 
3 m.-nouj 3 uj 


The above considerations give us the schematic time history of 
the electrostatic field energy as drawn in Fig. 11. Figure 11 
indicates a triple-peak structure in the electrostatic field 
energy. When the plasma ^ is small, 33 , a.^i are close. In this 
case, the triple-peak structure in the electrostatic field energy 
would become double-peak structure. The maximum value of the 
electrostatic field, Ep*v« , achieved at a = a^t is given by 

1 . 3 . ^ m- 

E. 


- 1 n 

a c® A ~ a X’ 


The induced electric field Ez 
that E becomes zero, when a 
a = 0 because a = 0 is the 

achieves maximum or minimum. 


is given by Eq. (30), which 
= 0 . Ei changes its sign 
point where the magnetic 


(45) 

shows 

around 

field 


Next, the time behavior of ion temperature T;x is examined. In 
the early stage of coalescence, the plasma should be 
ad i aba t i cal ly compressed. However, as the magnetic field energy 
increases near the peak and approaches the peak, the ion flow 
energy becomes dominant over the thermal energy. From the 
consideration that v^ gives maximum or minimum, namely^^ = 0 , 
we find two conditions for the extrema; (i) v^ = 0 , which gives 
a, (ii) = 0 , which gives 


aa = 


(46) 


When the explosive 
the condition (46) 


coalescence takes place 

„ 2 cl 

3 vr 


( £ = 0 ) , 


we estimate 
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1 1 . The temporal behavior of the electrostatic field energy constructed from the 
Sagdeev potential. 



12. The temporal behavior of the ion-temperature in the x-direction (in the 
direction of coalescence) constructed from the Sagdeev potential. 



After a = a^, the kinetic energy must decrease, which means that 
the plasma is in the state of colliding phase. The 
above considerations give us the schematic time history of the 
ion temperature, which is shown in Fig. 12 . Figure 12 shows a 
double-peak structure in the ion temperature. In the limit of 
quas i -neu t r a 1 i ty , we can estimate the dominant term governing the 
explosive phase where the adiabatic compression is predominant. 
The temperature T is given by T = P / n , while the dominant 
term in pressure changes in time as a"^ when V = 3, P — a"*^ when 
r = 2, while n a“* . Therefore we find 

T = P/n ^(-7 = 3) 

c. 1/a^^ = 2). 


(47) 

(48) 


We investigate in more detail the explosive phase of the 
coalescence in a case that we can neglect the effect of plasma 
pressure: it only acts as a saturation mechanism. However, if we 
take ^ = 2, the pressure does not take the role of a saturation 
mechanism as seen in Eq. (27). In the explosive phase, therefore, 
we can neglect the second term in the right-hand side of Eq. (26); 


Fu r the rmore , 
namely the 
examp 1 e , if 
order of the 
total energy 


a = 


v: 


(49) 


we need a second condition for explosive collapse, 
initial total energy must be close to zero. For 
the oscillation period of magnetic energy is of the 


transit Alfven time t/\ 

£ by making use of Eq. (40) 


we can estimate the 


as 


~ tA, 


(50) 


which gives 



(51) 


The solution of Eq. (49) with small £ is given by 

a(0 £ - f)2/3 + 0(f), (52) 

where we neglect the order of £ , tj, is the explosion time. Once 
the solution a (t) is given by Eq. (52), we can find the various 
physical quantities as follows, which is valid in the explosive 
phase of the coalescence; 

2 X 


n = n, = Tig 


/ 2 . 1/3 


(54) 
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(55) 


E^ = 


2 m, X 


,2.2/3 BooA‘/^x 

E _2/2v2/3 BqqA^/^x^ I . 

' 3^9^ _ f)7/3 3^9^ o;|,A»/3va/^to-0®/^’ 


Boo<^ 


initial 
at t = 0 


The explosion time t<, is related to the 
From Eq. (53) we find the initial velocity 

Ao 2x 

Vzo = — ^ 

Aq «JtO 

a„ are the initial values of a, a at t 


where a^ 

Eq. (58) we find 




(56) 

(57) 

condi t ion. 
as 


(58) 
0. From 


(59) 


where a- must be negative when magnetic collapse occurs. On the 


other hand, a and a are related to the 


€ = il- 


2v\ 


when the pressure term is neglected. 
Eqs. (59) and (GO) 

v/2 rf 

to — — flOV^O^A- 

O 


initial total energy as 

(GO) 


If £— > 0, we obtain from 


(Gl) 


Combining eqs. (55) and (5G) , we find an important result that in 
the explosive phase the electrostatic field Ex^(te “t) more 
rapidly grows compared with the magnetic fiel d (t^ -“t) . 
This fact becomes very important when we consider the high energy 
particle production by E during the explosive phase of the 
coalescence. 

Let us comp are the theoretical results obtained here for the 
explosive phase with the computer simulation results. The global 
structure of the magnetic field energy, electostatic field 
energy, ion temperature in the x~direction observed in the 
simulation is well explained by the theoretical model obtained 
here. Especially, the doubl e—pe ak structure in the ion 
temperature, and the triple -peak structure in the electrostatic 
field energy are also observed in the simulation (see Figs. 1 and 
2). Table 1 summarizes the results of comparison between the 
theory and the collisionless simulation. 

In Table 1 we show the index m of explosiveness (the 
exponent to the time (t^ -t)*”*^ ). Table 1 shows a good agreement 
between simulation and theory in the electrostatic energy. When 
the ttoroidal field becomes comparable the poloidal field (B-t*/ Bp 
= 1 case), there appears two-dimensional motion, plasma 
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rotational motion , which makes a more complex electrostatic 
field configuration. When Bt/ Bp = 1 , the magnetic 
time behavior of field energy also deviates from the one- 
dimensional analysis for the same reason. For the ion temperature 
we find its explosiveness from Eq. (47) as 




(G2) 


when 


T - 3. The above scaling also is close to the results 

obtained in the simulation, except for the case of Br/ Bf = 1. 

Figure 13 shows the spectral intensity SCk,(U ) obtained from 
the simulation run that corresponds to Fig. 2. Since the spectral 
intensity is integrated over time, the most intense period of 
various modes is most heavily weighted. For most of modes this 
intense period corresponds to that after the coalescence. Thus, 
it becomes necessary to recalibrate the plasma density and 
magnetic field strength, etc. at that stage. These recalibrated 
e i g e n f r e qu e nc i e s are indicated near the edge of the frame of 
Fig. 13 . We observe that the lower hybrid range frequency 

contains strongest spectral intensity, followed by 
c 0 mp r e s s i 0 na 1 Alfven modes that correspond to the coalescence 
ensuing oscillations, and the^plasma oscillations 
frequency happens to be outside of the frame so that we do 
know if they are prominent or not. 

Finally, comments in the effect of no n~qua s i neu t r a 1 i t y 

phase, we can neglect the effect 
and Eq. (24) : 


due. In 
pressure 


the 

terms 


explosive 
in Eq. (23) 


the 
and 
The 
no t 

are 
0 f 





,2 ’ 
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(63) 


(64) 


From the analysis of numerical calculations of Eqs. (63) and (64), 
we can conclude that b is slowly varying in time during the 
variation of a. Therefore if we use the result of b — be>-const. in 
Eq. (63), we find the effective potential V (a) as 

dV(a) 


a = —■ 


V(a) 


da ' 

= uz. I — a I — 


Me 


(65) 


14. 
0 f 


^‘[2bo V .A^a’ 
the potential is given in Fig. 
potential including the effect 
that incorporates the pressure effect is 
line near a = 0. Figure 14 shows a second 
is caused by the effect of charge 


Here we 
p 1 a sma 


The schematic graph of 
write the effective 
pressure. The curve 
shown as the broken 
minimum at a 

separation. When the J x B force dominates, the charge separation 
effect is unimportant. Such is the case of explosive phase of 
coalescence (where charge neutrality is maintained to a good 
degree). However, after the coalescence in the late stage of 


ap,!n which 
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amplitude oscillations, the J x B force becomes weaker. Then the 
charge separation may become important. In this stage electrons 
can be trapped and oscillate in the potential well near the 
second minimum in Fig. 14 . The electrostatic field oscillations 
are similar to the dipole oscillations near the current sheet. 
Such dipole oscillations may be able to radiate the 
electromagnetic wave C or ) , if the plasma is not 

evanescent. The simulation results obtained show certain high 
frequency oscillations in the frequency range of , Ct>Ln . The 

detailed compariso.n will be reported elsewhere. Once other 
effects such as the two-dimensional curvature effect are included 
, it might be possible to have additional minima in the effective 
potential and thus for the system to be temporarily trapped in 
the potential well and exhibit pulsations. 

We have found that in the explosive phase of magnetic collapse 
the electrostatic field can be explosively generated and more 
rapidly grow compared with the magnetic field. 

The explosiveness of the electrostatic field 
Ey (tp - t)*^ , and magnetic field 'v (t^- well agree 

with the results obtained in the simulations (see Fig. 1 ) 

As shown later, the kinetic simulation finds that in the 
explosive phase ions and electrons are simultaneously accelerated 
in the z-direction, opposite each other. When a particle moves 
across the magnetic field driven by the electrostatic field E 
the particle can be accelerated in the direction (z-direction) 
perpendicular both to the electric field (x-di rec 1 1 on) and the 
poloidal magnetic field (y-di rect ion) . This acceleration 
mechanism was considered by Sugihara and Mizuno (1979). On the 
other hand, Sagdeev and Shapiro discussed the same physics in 
another point of view, namely the large amplitude wave damping 
due to trapped electrons. These previous works are applicable to 
the cases in which the large amplitude electrostatic waves can 
propagate across the static magnetic field. In the explosive 
phase during magnetic collapse, a similar situation can be 
realized; now the magnetic field and electrostatic field can vary 
in time and space. 

Figure 15 shows the distribution functions of electrons and 
ions in the coalescence simulation from the electromagnetic 
particle code. In Fig. 15 we compare the distribution of 
electrons and ions before the explosive coalescence with the 
distribution after that. 

It is clear that the distribution functions have very rapidly 
changed during the explosive coalescence and strong heating of 
ions in the coalescence direction has occurred during this time. 
It is also observed that there is a very small but energetic 
population of electrons has been created in the z-direction 
perhaps due to the inductive acceleration. A similar pattern is 
observed in the ion distribution in the z-direction with the 
appropriate parity difference because of the charge difference 
between electrons and ions. It should be noted that the explosive 
coalescence produces extremely energetic ions which form a 
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relatively flat and long plateau. 


3. Current Loop Coalescence In Solar Flares 

Recent observations of x-ray continuum emission, T —ray line, 
and continuum emission from solar flares with instruments on the 
Solar Maximum Mission (SMM) and Hinotori satellites show that 
energetic ions and relativistic electrons are accelerated almost 
simultaneously with non-re lat ivi st ic electrons during the 
impulsive phase of solar flares. These observational results make 
it necessary to revise the widely accepted hypothesis of particle 
acceleration that energetic ions and relativistic electrons are 
produced in the second phase a few minutes after the impulsive 
phase (Wild et al. , 19G3; de Jager, 1969; Svestka, 1976). 
Although Bai and Ramaty (1979), Bai (1982), and Bai et al. (1983) 
revised the hypothesis as the secondstep acceleration taking note 
of a small delay of T -ray line emission from hard x-ray 
emission, Kane et al. (1983), and Forrest and Chupp (1983) 
pointed out that such a small delay can be explained simply by 
either the injection, propagation, or energy loss processes of 
particles which are accelerated in a single step. 

Recently Nakajima et al. (1983) and Kiplinger et al. (1983) 
reported observations of quasi-peri odic pulses with double sub- 
peak structure seen in hard x-ray, T -ray and microwave emissions 
in the two intense solar flares of June 7, 1980 and June 21, 
1980. We are interested in the close similarity between the 
observed time profiles and those obtained with the computer 
simulation given in the previous section. 

We present the results of our analysis of the June 7, 1980 and 
November 26, 1982 events, both of which are widely different from 
each other in duration, source size, source height, etc., they 
provide a stringent test for examining the validity of our model 
of particle acceleration in solar flares in terms of the 
coalescence instability. Our study shows that 
observational features of the two events are consistent with the 
results of our computer simulation. 


3. 1 Summary of Observations 
(a) June 7, 1980 Event 

The impulsive burst of the June 7, 1980 solar flare (Fig. 16) 
has been investigated by many authors (Forrest et al. , 1981; Kane 
et al. , 1983; Forrest and Chupp, 1983; Nakajima et al. , 1983; 
Kiplinger et al., 1983). We summarize below some essential points 
from these observations. 

(1) The burst is composed of seven successive pulses with a 
qua s i — pe r i o d i c i t y of about 8 seconds. Each of the pulses in hard 
x-rays, prompt -ray lines, and microwave is almost synchronous 
and similar in shape. 

(2) The microwave pulses consist of double sub-peaks as seen 
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1980 June 7 



Universal Time 


16. Time profiles of 17 GHz microwave and 150-260 kev X-ray emissions (from 
the Hard X-ray Burst Spectrometer on SMM) of June 7, 1980 Flare. 



17. (a) time profiles of 17 GHz microwave and (b) time variation of the height of 
two microwave (17 GHz) sources during the 1982 November 26 Flare, 
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especially in the second and fourth pulses in Fig. 16(a). 

(3) The starting times of hard x-rays, prompt T -ray lines, 
and microwaves coincide within ±2.2 seconds. 

(4) The time scales of acceleration for both electrons (up to 
energies above 1 MeV) and ions (above 10 MeV/nucleon) are less 
than 4 seconds. The accelerations must occur almost 
simultaneously. 

(5) The height of the microwave source is estimated to within 
10 are sec above the photosphere (H,< flare; N12 , W64 ). The 
source has a small size of less than 5 arc sec in the east-west 
direction and shows no motion. 

The Ho( photographs from the Peking Observatory (H. Chow, 
private communication) add a new finding. The flaring region has 
two structures that appear to be in contact with each other, one 
stretching in the east-west direction and the other in the north- 
south. 

(b) November 26, 1982 Event 

Ve briefly outline the characteristics of the Nobember 26, 
1982 flare (Fig. 17 ). This event is of much longer duration than 
the event on June 7, 1980, about 20 minutes compared to about 1 
minute. The microwave observations were made with the 17-GHz 
interferometer at Nobeyama, Japan, and the hard x-ray observation 
with the hard x-ray Burst Spectrometer (HSRBS) on SMM. 

(1) The microwave burst is composed of three successive peaks 
with a quasi-periodicity of about 6 minutes as indicated by 
number 1-3 in Fig. 17(a). 

(2) Each of the microwave peaks further consists of two sub- 
peaks. The hard x-ray time profiles seem to coincide with the 
microwave sub-peaks. The SMM hard x-ray data are available only 
for the first peak. 

(3) The microwave and hard x-ray emissions start almost 
simultaneously within 10 seconds. 

(4) The microwave source is composed of two sources, one at a 
height of ■^10^ km and the other at 3x10^ km. These values are 
derived on the assumption that the sources are located directly 
above the flare (SIO, W88) . 

Figure 17(b) shows the height of the two microwave sources as 
a function of time. In the pre-burst phase (phase 1: 0220-0228 
UT) , the upper source appears at a height of 2.9x10^ km above 
the photosphere and the lower one at 0. 7x10^ km. In phase 2, 
the lower source rises at a velocity of 30 km/ s . The main 
phase (phase 3) started when the lower source reaches a height of 
1. SxlO^km. It is suggested that the two sources collide with each 
other at this time. In fact, a small up-and-down motion of 
the lower source is observed in the main phase. The oscillation 
period and peak-to-peak amplitude of the up-and-down motion are 
/''min and <'^2x10^ km (significantly larger than the fluctuation 
level due to the signal to noise ratio), respectively. After the 
main phase, the lower source begins to go down towards its 
previous position. On the other hand, the upper source rises 
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gradually, though it remains at almost the same height until the 
decay phase. 

The observational facts summarized above, especially the 
collision of the two microwave sources and the small up“and~down 
motion of the lower source in the November 26,1982 event, suggest 
that the current loop coalescence takes place. The existence of 
two bright components in the June 7,1980 event also supported 
this interpretation. 


3*2 Interpretation by Simulations 

Two parallel current loops are unstable against the coalescence 
Instability. They are attracted by and collide with each other 
and finally coalesce into one loop. Its nonlinear development can 
release a large amount of poloidal magnetic energy associated 
with the current loops into particle energies. We investigated 
this process, i.e.,the global plasma dynamics, heating and 
acceleration of pa r t i c 1 e s , and so on, through computer simulations 
as described in Chapter 2. Here, we made two different types of 
simulations in order to experiment with a wide variety of plasma 
parameters: one is a MHD particle simulation, and the other a 
collisionless full-electromagnetic particle simulation, both of 
which are two-dimensional in space across the plane perpendicular 
to the current loops and three dimensional in velocity space. 

(a) Explosive Coalescence June 7,1980 Flare 

The case that two parallel loops have sufficient electric 
currents so that they attract each other fast enoughCin about one 
Alfven transit time) was simulated using the collisionless full- 
electromagnetic particle code. 

The resultant time history of the electron temperature is 
shown in Fig. 2. We can clearly see a qu a s i -p e r i o d i c oscillation, 
the period of which is about one Alfven transit t i me (8 ) . The 
cause of this oscillation is as follows: after explosive 
reconnection of poloidal magnetic fields takes place at the X- 
point between approaching current loops, the two plasma blobs 
pass through each other and overshoot, resulting in repetition of 
this process. 

Fig. 2 also shows that the electron temperature oscillation is 
characterized by prominent double sub-peak structure. The double 
sub -peaks occur just before and after each peak in the magnetic 
field intensity. Just before a peak, the magnetic acceleration of 
the plasma by JxB becomes strongest so that the magnetic flux 
behind the colliding plasma blobs as well as the plasma blobs 
themselves are strongly compressed. This plasma compression 
causes the first sub-peak of the electron temperature. Then, the 
plasma particles acquire velocities close to the Alfven velocity 
along the colliding direction, so that they detach from the 
magnetic flux against which they have been compressed, resulting 
in an expansion and hence in an adiabatic cooling of the plasma 
as the magnetic fields obtain peak values. After the peak in the 
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magnetic fields, the process reverses giving rise to the second 
sub-peak of the electron temperature. 

A similar time history is obtained for the kinetic energy of 
high-energy tail electrons and protons as well as for proton 
temperature. The acceleration of the high energy-tail particles 
is due to the Vp xB acceleration mechanism, which can lead 
simultaneous electron and proton acceleration within about 
second. Since these processes accompany the local plasma 
compression/decompression just before and after coalescence, it 
is not surprising that the time profile of the microwave 
emissions caused by high-energy tail electrons resembles that of 
Fig. 2(d). 

The result of this simulation can also explain the observed 
period of the quas i-pe r i od i c oscillation of the June 7, 1980 

event. The period( 4 sec) which is estimated with source size 
arc sec), magnetic fieldC'^200 Gauss: Kiplinger et al.,1983) and 

emission measure (/%. 10*"^ cm"'^ f rom the GOES soft x-ray data). 

We discuss the energy spectrum of electrons and protons, after 
the explosive coalescence of two current loops. The energy 
spectrum consists of three components ; (a) background thermal 
component due to the adiabatic heating, (b) i nt e rmed i a t e component 
due to inductive electric field, (c)high energy component due to 
the Vo xB acceleration. Fig. 18 shows the energy spectrum without 
the background thermal component which was reproduced from the 
previous figures given in Chapter 2. 

The intermediate non-thermal component has V- 2 as the power-law 
index near the peak, while near the valley the spectrum becomes 
more soft. The global structure of the electron energy spectrum 
is consistent with observations. (Kane et al. , 1984; Kiplinger et 

al. , 1984) 

The spectrum for ions is shown in Fig. 18b in which two cases are 
given. The one is just after the coalescence( T=28ils,*) and the 
other is at T = 40J1;' . As seen in the figure, the spectrum 
becomes harder which means that the number of high energy proton 
increases. The range (b) and (c) for protons corresponds to the 
production of observed T -rays. 

Finally, we make brief comments about the energy spectrum in 
the multiple coalescence. The global behavior consisting of three 
energy components is still the same as the two loop coalescence 
(see Fig. 19 ). The maximum energy obtaind by Vp x B acceleration 
depends on total released energy. Therefore the bigger flare can 
produce stronger -rays and neutrons. 

(b) Slow Coalescence November 26 1982 Flare 

When two parallel loops have insufficient electric currents or 
are not well separated and hence the attracting force of them is 
weaker than that of the previous case, reconnection of poloidal 
magnetic fields during loop coalescence becomes slower. (However, 
this reconnection rate is still faster than what would be 
predicted by a classical tearing theory (Furth et al. , 1963)). 
This case was also simulated using the MHD particle code. Figure 
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18. (a) Electron energy spectrum just after explosive coalescence, (b) Proton 
energy spectrum just after (t=28) and later explosive Coalescence (t=40). 
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19. Distribution functions for electrons and protons in the case of multiple-current 
loop interactions. 
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7 (a) shows the temporal development of plasma kinetic energy 
Cthe electron pressure energy) during the coalescence. Also shown 
in Fig. 7 (b) is the time history of the integrated reconnected 
magnetic flux through the X-point (case shown in Fig. 5(a) in 
Bhat tachar j ee , et al. , 1983) Note that a slight amount of 
oscillations of reconnected flux can be seen around the straight 
line. Again we can see the oscillatory behavior with double sub- 
peak structure in both time histories, though it is less 
prominent compared with that of the explosive coalescence case 
presented in the previous subsection. The period of oscillation 
is about 5 times the Alfven transit time. 

The obtained time history resulting from the simulation is 
explained as follows. In the case of slower reconnection, the two 
plasma blobs do not pass through each other but are pushed back 
by the pressure of the magnetic field compressed between the two 
loops. This motion is repeated resulting in the damping 
oscillation shown in Fig. 7 (a). The amplitude of the oscillation 
in this case is less prominent compared with the previous case. 

The observed plasma kinetic energy oscillation exhibit a 
structure quite similar to the microwave time profile of the 
November 26, 1982 flare as shown in Fig. 17(a). The source size of 
the November 26 flare is about 10 times larger than that of the 
June 7 flare. We therefore estimate the calculated period of the 
oscillation to be 5 x 4 x 10 = 200 sec, assuming that the Alfven 
velocity is about the same for both cases. This period is close 
to the observed period of about 6 minutes. Note also that in this 
case the flow velocity is much below the Alfven velocity in 
agreement with the observational fact that the 40 km/s colliding 
velocity of the lower loop is much smaller than the Alfven 
velocity of 10^ km/s. 

The results obtained from computer simulations of the 
coalescence instability of two current loops are in good 
agreement with observations of two widely differing flares. The 
key characteristics which are well explained are the simultaneous 
accelerations of both electrons and ions, and the double sub-peak 
structure in quasi-periodic pulses. The double sub-peak structure 
is more pronounced when the currents in the two loops are 
sufficient for the explosive coalescence to occur. This case 
corresponds to the June 7, 1980 flare. When the currents are 
insufficient for the fast coalescence, the double sub-peak 
structure is less pronounced. This case corresponds to the 
November 26, 1982 flare. In addition, we have the observation 
suggesting the collision of the two microwave sources for the 
November 26,1982 event. 


429 


4. Conclusions and Discussion 

In the preceding chapters we have shown fundamental 
characteristics of current loop coalescence plasma dynamics and 
compared them with solar flare observations. The current loop 
coalescence has been successful in explaining several important 
natures of the impulsive phase as f o 1 1 ows : (1) the sudden explosive 
development in time p r o f i 1 e s , (2) s i mu 1 1 ane ou s acceleration of 
electrons and protons up to several times their rest mass energy, 
C3)plasma heating due to adiabatic compression up to several 
times 10 K, and (4) quas i ~pe r i o d i c pulses of microwave, x~ray and 
f-ray radiations and also double-peak structure in each pulse. 
With regard to C2) , it is noteworthy that the enegy spectrum of 
accelerated electrons is almost consistent with observations. 
Further, we have revealed that the energy release due to the 
coalescence can be very explosive or relatively gradual depending 
upon various physical conditions, the key parameters of which are 
plasma ratio, colliding velocity (current localization and 
peakness condition), and Bp/ • Thus this model is applicable 
both to very ’’impulsive flares” and to ’’gradual flares”, as is 
shown in Chapter 3. Probably in impulsive flares such as the June 
7, 1980 flare small current loops develop and interact with each 
other in the lower corona, while in gradual flares such as 
November 26, 1982 flare, large current loops whose to reaches 
several 10 km in its height, develop and interact with each 
other, whose current distribution might be comparable with the 
distance of the two loops. 

So far, we have argued that the quas i -pe r i od i c i t y is one of the 
most important signatures of the two current loop coalescence, 
because of its ideal situation. However, the qu a s i -p e r i o d i c 
pulsation is not a common phenomenon for the majority of flares. 
In particular, the beautiful regularity such as seen in the June 
7, 1980 flare has been rarely observed. Therefore one might claim 
that the current loop coalescence in not a basic mechanism of 
solar flares, even though in some specific flares it plays a 
role. On the contrary, we believe that the coalescence is an 
elementary process commonly occurring in the impulsive phase of 
solar flares. 

In Chapter 2 we have presented the result of a computational 
simulation of multiple interaction of many current filaments in 
addition to that of two current loop interaction. The multiple 
coalescence exhibits many bursts that are more irregular than the 
simple two current loop coalescence. In other words, it is very 
difficult for us to find any regularity from the resultant time 
profiles alone. We suggest that this is the case that happens in 
most flares. As is easily expected, two dominant loops with 
comparable current with each other may be (rarely) hardly 
realized in f 1 a r e -p r o due i ng active regions where the photospheric 
magnetic configuration and the photospheric motion are very 
comp 1 i cated. 

As shown in the previous sections, the current loop 
coalescence model has been successful in various observed points 
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in the impulsive energy release stage. We have shown the many 
characteristics of the coalescence dynamics as the element 
process of current -current interactions which might be basically 
of fundamental importance for the energy release process. We now 
discuss the preflare stage in which multiple current-filament 
structure as shown in Fig. 20 might be generated due to the 
photospheric shear motions. The observations near the active 
region imply that before several hours of the onset of the flare 
there exists photospheric shear motion across the magnetic 
neutral line as shown in Fig. 20 . The photospheric shear motion 
can give rise to plasma currents along the potential magnetic 
field produced by the sunspots near the active region. If the 
photospheric motions are associated with plasma vortex motion, 
the current loop structure can be generated in the arcade-like 
structure as shown in the figure, 

As the shear motions proceed, the current density can increase 
and the current loops might move up> associated with relaxation 
of magnetic tension. Recent computer simulation (Wu, 1985) 
suggests in the modeling of arcade-like preflare stage that the 
current tends to localize by its pinching effect in the low 
region. Recently, the preflare acceleration in the current sheet 
moving up across the external potential magnetic field is 
discussed (Sakai et al., 1985). Furthermore, the recent 
simulation (Aydemir et al., 1985) of current generation in the 
current loop with line-typing effect suggests the multi-current 
filament structure with different helicity as shown in Fig. 21. 
The multi-current filament system might ascend with increase of 
the current density as well as current constriction into the low 
region of the corona. In the low (i region the current 
constriction can be enhanced and each current filament can be 
well separated. If the current density continues to increase 
further, such that the condition Bp >Bt can be realized, the 
current loop coalescence can be set up. Of course, as se.en in the 
mu 1 1 i -c o a 1 e s c e nc e simulation, the not-well separated current 
filaments might coalesce without a large amount of energy release 
and can grow to high current density filaments before explosive 
strong energy release, A similar situation can also occur in the 
braized current filament system in Fig. 21. 

Next we discuss the later phase after the current loop 
coalescence in connection with the observed two Ribbon structure. 
As discussed above, the flaring region might consist of multiple 
current loops associated with ascending motions. As shown in the 
figure, if the current loop coalescence happens to start locally 
where the condition for the explosive onset can be satisfied in 
the lower current loops, the sequential flaring might precede 
firstly to the direction along the magnetic neutral line, as 
observed in the recent Hot observation (Kurokawa et al. , 1985). 
The reason why the sequential flaring proceeds along the arcade 
might be due to geometrical structure or the fact that the meta- 
stable filaments sit around in the same height in the corona. The 
flaring loops move up slowly and continue to trigger the current 
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loop coalescence with upper currect loops. The sequential flaring 
on the upper loops can be observed in the picture which shows 
slow expansion motions far from the magnetic neutral line. 

We need more investigations about current filament generation 
in the preflare stage, in connection with the explosive current 
loop coalescence process. 
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